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(57) Abstract 

An apparatus and method for computing both tomographic and panoram- 
ic scans of an object (16). The apparatus (10) includes an X-ray source (12) and a 
linear x-ray detector (20), and is operable in two modes: a panoramic mode and 
a tomographic mode. In the panoramic mode, the detector (20) is aligned per- 
pendicular to a plane of rotation, and the detector (20) and x-ray source are ro- 
tated together about an object (16) to be imaged, and about a center of rotation. 
In addition, the center of rotation travels along a predetermined path. The digital 
signals produced by the detector (20) are electronically processed to produce a 
panoramic x-ray image of the object (16). In the tomographic mode, the linear 
detector (20) is positioned to lie within the plane of rotation, and the detector 
(20) and x-ray source (12) are rotated together in the plane of rotation about a 
center of rotation and about an object (16) to be imaged. Signals produced by the 
detector (20) during rotation are electronically processed to produce a computed 
tomographic image of the object (16) being scanned. During the tomographic 
mode, the center of rotation can also travel along a predetermined path, and this 
eccentricity is mathematically corrected. 
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METHOD AND APPARATUS FOR COMPUTING 
TOMOGRAPHIC AND PANORAMIC SCANS 



The present invention relates to an apparatus and 
method for producing computed tomographic (CT) x-ray 
images as well as panoramic x-ray images. 

X-ray computed tomography (CT) is a technique for 
obtaining cross-sectional reconstructions of three 
dimensional objects using x-rays. In the simplest 
example of CT imaging, a narrow beam of penetrating 
x-rays is scanned across an object or patient in 
synchrony with a radiation detector on the opposite side 
of the patient. If the beam is monoenergetic or nearly 
so, the transmission of x-rays through the patient is 
given by the equation 

I = I 0 exp(-Mx) C 1 ] 

where the patient is assumed to be a homogeneous medium 
with the attenuation coefficient p. If the x-ray beam is 
intercepted by two regions with attenuation coefficients 
M, and n 2 and thicknesses x x and x 2 , the x-ray transmission 
is characterized as 



I = I 0 exp - (MiX t + M2X2) 



[2] 



This formula is generalized to many (n) regions with 
different linear attenuation coefficients with the 
argument of the exponent 



i=l 



1 = EI„exp - (m) t4] 
i=l 



separate attenuation coefficients cannot be determined 
with a single transmission measurement because there are 
too many unknown values of ft in the equation. However, 
with multiple transmission measurements at different 
orientations of the x-ray source and detector, the 
separate coefficients can be distinguished so that a 
cross-sectional display of coefficients is obtained 
across the plane of transmission measurements. By 
assigning gray levels (see below, to different ranges of 
attenuation coefficients, a display is obtained that 
represents various structures in the patient with 
different x-ray attenuation characteristics. This gray 
scale display of attenuation coefficients constitutes a 
CT image. 

The numbers computed by the reconstruction algorithm 
are not exact values of attenuation coefficients 
instead, they are integrals, termed CT numbers, which are 
related to attenuation coefficients. On most newer « 
units, the CT numbers range from -1,000 for air to *1.000 
for bone, with the CT number for water set at 0. CT 
nral bers normalized in this manner are termed Hounsfield 
units and provide a range of several CT numbers for a one 
percent (1%) change in attenuation coefficients. 
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To portray the CT numbers as a gray scale visual 
display, a CRT is used. The CRT includes a contrast 
enhancement feature that superimposes the shades of gray 
available in the display device (i.e., the dynamic range 
5 of the display) over the range of CT numbers of 

diagnostic interest. Control of image contrast with the 
contrast enhancement feature is essential in x-ray 
computed tomography because the electron density, and 
therefore the x-ray attenuation, is remarkably similar 
10 for most tissues of diagnostic interest. These electron 
densities vary from 3.07 x 10* elec/cc for fat tissue to 
5.59 x 10 23 elec/cc for the densest tissue, bone. Lung 
tissue has a much lower electron density, 0 . 83 x 10 23 
elec/cc, because of the alveolar and bronchial spaces. 

15 

The first CT systems were introduced in 
approximately 1971 by the EMI Corporation of England. 
These early systems used an x-ray source mounted in a 
gantry with detectors. The patient was inserted between 
20 the x-ray source and the detectors and the joined x-ray 
source and detectors were moved about the patient to 
obtain projection rays through the patient. These values 
were fed to a computer which then reconstructed a cross 
sectional image of the plane through which the pencil 
25 beam of x-rays passed. During this translational scan of 
perhaps 40 cm in length, multiple (e.g., 160) 
measurements of the x-ray transmission were obtained. 
Next, the angular orientation of the scanning device was 
incremented one degree and a second translational scan of 

30 160 transmission measurements was performed. This 

process was repeated at one degree increments through an 
arc of 180 degrees so that 28,800 x-ray transmission 
measurements were accumulated. Those measurements were 
then transmitted to a computer equipped with a 

35 mathematical algorithms for reconstructing an image of 



attenuation coefficients across the anatomical plane 
defined by the scanning x-ray beam. 

Although this approach yielded satisfactory images 
of stationary objects, considerable time (4-5 minutes) 
was required for data accumulation and the images were 
subjecTto motion blurring. Soon after the introduction 
cf pencil beam scanners, fan-shaped x-ray beams were 
introduced so that multiple measurements of x-ray 
transmission could be made simultaneously. Fan beam 
geometries, with increments of a few degrees for the 
different orientations (..,.. a 30-degree fan beam and 
10-degree angular increments) , reduced the scan time to 
20 -60^econds and improved the image quality by reducing 
the effects of motion. Computed tomographic scanners 
with x-ray fan beam geometries and multiple radiation 
detectors constituted the second generation of CT 
scanners . 

in late 1975, the third generation of CT scanner was 
introduced. These scanners eliminated the translation^ 
motion of previous scanners, using rotational motion of 
the x-ray tube and detector array or rotational motion of 
Z x-ray tube within a stationary circular array of 600 
or more detectors. With these scanners, data 
accumulation times as fast as two seconds are achievable. 

Both stationary and rotating anode x-ray tubes are 
used in CT scanners. Hany of the translation— rotation CT 
scanners have an oil-cooled, stationary anode x-ray^ube 
with a focal spot on the order of 2 X 16 mm. The limited 
cutput of these x-ray tubes necessitates a sampling time 
cf about 5 msec for each measurement of 

transmission. This sampling time, together with the time 
required to move and rotate the source and detector. 



limits the speed with which data can be accumulated with 
CT units using translational and rotational motion. 



To reduce the sampling time of 2-3 msec, most fast- 
scan CT units use rotating-anode x-ray tubes, often with 
a pulsed x-ray beam, to achieve higher x-ray outputs. 
Even with rotating-anode tubes, the heat-storage capacity 
of the anode may be exceeded if cooling periods are not 
observed between sets of successive images. 

After transmission through the patient, the x-ray 
beam is collimated to confine the transmission to a slice 
with a thickness of a few millimeters and to reduce 
scattered radiation to less than one percent (1%) of the 
primary beam intensity. The height of the collimator 
defines the thickness of the CT slice. This height, when 
combined with the area of a single picture element 
(pixel) in the display, defines the three-dimensional 
volume element (voxel) in the patient corresponding to 
the two-dimensional pixel of the display. A voxel 
encompassing a boundary between two tissue structures 
(e.g., muscle and bone) yields an attenuation coefficient 
of the two structures. This "partial volume artifact" 
may be reduced by narrowing the collimator to yield 
thinner slices. However, this approach reduces the 
intensity of the x-rays incident upon the detector and 
the detector signals are subject to greater statistical 
fluctuations, thus introducing more noise in the 
displayed image. 

To reduce the detector response time, all detectors 
used in CT scanning are operated in current rather than 
pulse mode. Also, rejection of scattered radiation is 
assigned to the detector collimator rather than to pulse 
height analyzers. Detectors for CT scanning are chosen 
on the basis of detection efficiency (greater than 50%) , 



short response time and stability of operation, and are 
either gas-filled ionization chambers or solid 
scintillation detectors. With any detector, the 
stability of response from one transmission measurement 
to the next is essential for the production of 
artifactfree reconstruction images. With a pure 
rotational source and detector geometry, for sample 
detector instability gives rise to ring-shaped artifacts 
in the image. Minimum energy dependence of the detector 
over the energy range for the CT x-ray beam also is 
important if corrections for beam hardening are to be 
applicable to all patient sizes and configurations. 

All of the early CT systems were designed and built 
only to perform CT studies. The concept of using other 
types of radiation sources that had not been specifically 
Signed for CT imaging was initiated in the mid 1970's. 

several of these efforts utilized existing x-ray 
therapy simulators. An x-ray simulator is a device that 
duplicates a radiation treatment unit in terms of its 
geometric, mechanical and optical properties but uses 
diagnostic x-ray tube as the source of radiation to 
siaulate the properties of the treatment beam. A 
simulator allows the beam direction -d the treatment 
fields to be determined while encompassing the target 
object with the simulator's irradiation. Since the 
simulator's emissions are generally less intense and less 
energetic than the emissions of therapy devices thereis 
a reduction in the target object's exposure to radiation. 

The combination of a detector system and an x-ray 
therapy simulator provides the necessary front end of a 
CT system. Application of the requisite information 
processing techniques and algorithmic reconstruction 



processes in combination with the simulator cum detector 
system enable production of CT images. 

In conventional CT scanners, the x-ray source emits 
a scanned beam or fan of x-rays that lie in a plane, and 
the x-ray detector lies in the same plane. The x-ray 
source and detector are substantially fixed relative to 
one another, and are rotated together in the same plane 
that contains the x-ray detector and the x-ray beam. In 
contrast, in another type of radiography, known as 
panoramic radiography, the x-ray source and detector once 
again lie in the same plane, however, they are rotated in 
a plane that is transverse to the plane containing the 
source and detector. In certain panoramic x-ray devices, 
such as those used to produce panoramic images of the 
maxillo facial region, the center of rotation of the axis 
connecting the x-ray source and sensor moves along a 
curvilinear path while the sensor and detector are 
rotated about the center of rotation. During a panoramic 
scan, the detector, which is preferably a digital 
radiation detector, is periodically read, and the data is 
then conditioned to produce a digital panoramic x-ray 
image of the target object. In addition, and in contrast 
with conventional CT scanners, panoramic x-ray imaging 
devices typically scan only approximately 230<» to 240° 
about the target object. 

Reconstruction Algorithm 

The first reconstruction formula was derived by an 
Austrian mathematician, P.J. Radon, tjber die Bestimmung 
von Funktionen durch Ihre Integralwerte Laengs Gewisser 
Mannigfaltigkeiten, Berichte Sachsische Akadamie der 
Wissenschaf f en (Leipzig) , Mathematische-Physische Klasse 
69 262-277 (1917) . His methods were based on the 
assumption of use of a parallel beam and concomitant data 



collection. The fundamental process can be described as 
a convolution followed by a back projection. Subsequent 
methods derived from the Radon method utilized variations 
in the application of filters and convolution to correct 
for the effect of noise in the reconstruction process. 
Pavkovich and Nunan (U.S. Patents 4,149,248 and 
4 365,339) modified Radon's algorithm to accommodate fan 
beam geometry. Their method consisted of a convolution 
of the original projection data with a function derxved 
from the inverse Fourier transform (IFT) of a product of 
the Fourier transform of the projection data and the ramp 
function in frequency space. These results were then 
back projected to the point of interest, (X,Y) . Their 
convolution was essentially an angular convolution for 
any given point. All of their mathematical manipulations 
were performed in configuration (real) space with no 
resort to transforms to frequency space. 

Algorithms based upon the Radon method such as that 
developed by Pavkovich and Nunan suffer from two defects: 
(1) they are computationally burdensome, requiring both 
excessive memory and time of execution, and (2) they are 
very sensitive to noisy data; that is, resolution of the 
reconstructed image decreases disproportionately with 
increases in noise in the data. The Radon reconstruction 
method is theoretically perfect when operating on perfect 
data; that is, an infinite number of parallel, noiseless 
projections detected by an infinite number of detectors 
of infinitesimal dimension. These requirements 
unfortunately do not exist in physical systems. Using 
real data, results from the Radon method for calculating 
attenuation coefficients degenerate rapidly as the data 
becomes noisy and a finite number of discretely sized 
detectors is used. 



Known reconstruction algorithms are of four general 
types: 



(1) simple back projection (Radon); 

(2) integral equations; 

(3) Fourier transform; and 

(4) series expansion. 

In the present invention, the Fourier transform is 
utilized. In this approach, the x-ray attenuation 
pattern at each angular orientation is separated into 
frequency components of various amplitudes, similar to 
the way a musical note can be divided into relative 
contributions of different frequencies. From these 
frequency components, the entire image is assembled in 
••frequency space" and then reconstructed by an inverse 
Fourier transform reconstruction process into a spatially 
correct image. 

While the various algorithms used for CT image 
reconstruction each have their own limitations, the 
quality of the overall procedure is dominated by the 
quantity and quality of the measured transmission data. 
The quantity of data is restricted by the specific 
scanner design and by limitations placed on time and 
computer resources. Phenomena which tend to degrade the 
quality of the measured data include: 

a) Geometrical errors such as misalignment or 
motion of the scanning system or patient 
motion, 

b) Instability of the x-ray source, 



.-10- 



30 



c) 



Statistical fluctuation of the measured signal, 

d) Polychromaticity (non-monochromaticity) of the 
x-ray beam, 

5 e) The finite dimensions of the scanning aperture, 

f ) Residual signal due to the time response 

function of the detector system (afterglow) . 

10 if as a result of these factors, the projection values 
derived from the measured data do not adequately 
represent the line integrals of the linear attenuatxon 
coefficients within the slice being scanned, even the 

15 most perfect reconstruction algorithm will give rise to a 
distorted image. Each of these factors and the manner xn 
which they are addressed in the case of the present 
invention is set out in the following paragraphs. 

20 Geometrical Error 

The most practical approach to error introduced by 
scanner misalignment or motion lies in the construction 
and maintenance of the scanner itself, include periodic 
2S testing of the mechanical registration. Patient motion 
is less controllable but can usually be minimized by 
proper patient support and through the use of the fast 

'fan beam, scanner. If patient motion is -stored such 
as through the use of a transducer arrangement or a laser 
beam reflection method, then data correction is feasible. 
Using the method of the present invention, correction can 
even be made for eccentricity in the motion of the center 
of rotation; that is, «hen the center of rotation does 
not coincide with the center of reconstruction at a given 
3S range. Such eccentricity is measured and subsequently 
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corrected for all angles involved in the back projection 
and reconstruction. 

In addition, the present invention can accommodate 
intentional eccentricity in the motion of the center of 
rotation, for example, when a panoramic x-ray apparatus 
is modified, in accordance with the present invention, to 
conduct CT scanning. In such CT scanning the center of 
rotation can intentionally move along a curvilinear path. 
Even such a gross eccentricity can be measured and 
subsequently corrected for all angles involved in the 
back propagation and reconstruction, in accordance with 
the present invention. 

X-Ray Source Instability 

The instability of the x-ray source is generally 
corrected for by adjusting the measured data in 
accordance with the signal measured by a reference 
detector. Another approach is to monitor the electronic 
parameters of the x-ray source, such as kVp and mA, and 
to make the corrections to the measured data from this 
information. The reference detector method, while not 
extremely sensitive to kilovoltage variations, is simpler 
and more easily utilized through electronic hardware. 
Computer correction in accordance with the present 
invention makes possible the use of lower cost x-ray 
power supplies and kVp monitoring to correct for changes 
in effective beam energy. 

Statistical Fluctuation 

Systemic error may arise in such forms as drift and 
gain variations in the detector system and associated 
electronics, or in the form of background or bias 
currents. If these variations can be monitored and 



quantified, corrections can be implemented 
hardware circuitry or computer correction. 



Polychromatic Effects 

Since flux-rate requirements based on statistical 
considerations for any reasonable scan time normally rule 
out the use of other sources of radiation, an x-ray tube 
is the only practical photon source for CT scannxng. As 
a result, a spectral distribution of photons is involved 
rather than photons of a single energy. Since lower 
energy photons have higher attenuation coefficients than 
higher energy ones, the beam becomes progressively 
"harder" as it traverses an increasing patxent thxckness. 
A "harder" beam, having a lower effective attenuation 
coefficient than a "softer" one, introduces a degree of 
inconsistency into the measured data used for 
reconstruction. In the absence of compensatxon for thxs 
effect, e.g., utilizing a fixed length water bath or 
software corrections applied at the preprocessing stage, 
this effect leads to a distorted image characterized by a 
general increase in coefficients from the center to the 
periphery of the cross-section. 

The object of preprocessing the measured 
transmission data before reconstruction is to "linearxze 
the logarithm of the ratio of the incident to transmxtted 
intensity. For a monoenergetic beam traversing a 
homogeneous medium, this logarithm is a linear functxon 
of increasing thickness, while for a polychromatxc beam 
this function is no longer linear. If the chaxacterxstxc 
attenuation of a particular incident x-ray beam by an 
increasing mass thickness of water is known and xf xt can 
safely be assumed that the materials encountered wxthxn 
the body have attenuation properties similar to those of 
water, the measured data can be corrected to an xdealxzed 
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monoenergetic (linear) response for some suitable energy. 
Since the corrected data is then logarithmically linear, 
that data can be utilized by the reconstruction algorithm 
to produce a spatially consistent image that is 
5 independent of patient size. This correction may also be 
implemented by assuming an average composition of tissue 
and bone instead of water for the various degrees of 
attenuation. 

10 Finite Dimensions of Scanning Aperture 

The Qualities of resolution and image sharpness are 
closely associated with aperture size. Theoretical 
development of the reconstruction algorithms is based 

15 upon an infinite amount of infinitely thin transmission 
data which, in practice, is approximated by a finite 
number of transmission measurements of an x-ray beam of 
finite dimensions. The minimum aperture size is limited 
by photon statistics for a given x-ray output, geometry 

20 and scan speed. 

The slit height, perpendicular to the linear motion 
of the scanner, acts to make the resultant reconstructed 
coefficient a type of average coefficient over many thin 

25 transverse planes. This average is not inconsistent with 
the assumption of infinitely thin rays since the 
theoretical development of these algorithms is limited 
to two dimensions. Distortions due to this vertical 
smearing are somewhat minimized because of the 

30 homogeneity of the human body over short vertical 
distances. In some respects, this smearing is 
advantageous in that a larger volume is considered in 
each cross-sectional slice, hence fewer slices need be 
reconstructed to include the entire volume of interest. 

35 However, in the two-dimensional reconstruction, the 

resultant coefficient need not pertain to only one type 



of tissue, especially for small structures and near 
boundaries. This necessity is particularly important xn 
relating these coefficients to effective atomic number, 
density, chemical composition, or specific tumor types. 

The slit width, parallel to the linear motion of the 
scanner, introduces a compromise between the algorithm 
development and practical measurement. The approximation 
that the average intensity transmitted along the width of 
the aperture is the same as the relative intensity along 
a central ray is partially responsible for the lack of 
sharpness noted along many boundaries and can 
significantly degrade resolution. If the slit width is 
less than the linear increment between samples, computer 
enhancement of the data is not feasible. However, if the 
slit width exceeds this increment, then information is 
available from which the intensity for an aperture of a 
width approaching that of the linear increment can be 
calculated in accordance with the present invention. 
This technique is comparable to ^convolving the point 
spread function from other imaging devices, except the 
correction is applied to the measured data before an 
image is formed rather than as a modulation transfer 
function enhancement done as a post processing procedure 
on the final reconstructed image. The deconvolution is 
along discrete steps, the measured data, and not a 
continuum, therefore leading to the limit of data 
resolution, the linear increment. 

Time Response of Detector System 

The time response is an important consideration in 
the choice of detector systems. If signal decay do* to 
an impulse of radiation on a detector is slow, then it 
would appear from the measurements alone that some 
radiation was still incident on the detector some short 
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time after radiation exposure. Likewise, a particular 
measurement during a scan may be partially due to 
radiation incident on the detector from some time prior 
to that measurement. In the case of the linear scanner 
5 considered here, this temporal smearing due to the 

detector response may be related to a type of spatial 
response of the detector by the time spacing of the 
measurements or scan speed. In this way, the time 
response of the detector is considered similar to or as 
10 part of the aperture transmission function and is 
corrected for in the same manner. 

Applicability to Fan Beam Scanners 

15 For fan beam scanners, many of these correction 

techniques still apply. Furthermore, the temporal 
response of the detector is related to measurements made 
with the same detector as a function of angle (time) . 
The aperture itself is less likely to overlap than in the 

20 case of the linear scanner, however, the effective 

aperture due to cross-talk and patient scatter may also 
be handled by these preprocessing techniques. 
Modifications applicable to patient scatter may include 
making a correction to the measurement of interest as a 

25 function of intensity attenuated (including scatter) 

rather than intensity transmitted to a nearby detector 
along with the distance of this detector from the 
detector of interest. In any case, this type of scatter 
correction is just an approximation. The corrections 

30 pertaining to geometrical errors, instability of the 

x-ray source and the polychromaticity of the x-ray beam 
are applicable to either type of scanner. 

Some of these techniques are known to be in 
35 commercial use, especially in regard to polychromatic 

correction. Hardware approaches to some of the problems 



«. also common, such as careful construction of the 
scanner and use of patient supports and restraints to 
minimise geometrical error, reference detector methods to 
compensate for x-ray source variation, temperature 
compensating amplifiers to reduce drift, "f-^jT 
filtration to reduce polychromatic effect, ""-""V x 
ray output to allow for smaller apertures, and the 
selection of detectors to reduce afterglow Even so, the 
present method for reducing these factors which degrade 
L guality of the data is beneficial and in some cases, 
absolutely necessary. 

in spite of these many improvements and alternative 
approaches, there remains a need for improved methods for 
reducing the effect of these error factors on a back 
projected image, as well as for actually back P"3 acting 
L image of the target object. There is also a need for 
an apparatus for implementing these methods, and 
.specially an apparatus which is capable of being mounted 
in combination with an existing x-ray simulator Without 
changing the function of that apparatus. 

In addition, several fields of medicine exist which 
use both CT scanning and panoramic imaging, typically 
requiring two separate machines. There is thus a need 
fofa single x-ray apparatus capable of performing both 
CT scanning and panoramic imaging. 

.l.w ^ <Pir<-*- actiPct of the invention, 
In accordance with a first aspect or 

tt e present invention solves the above-noted drawbacks of 
the prior art by providing an attachment for back 
oroiecting a CT scan using an existing x-ray simulator 

. linear detector array, the mounting of which 
is adapted for mounting to the film holder of an x ray 
simulator normal to the central axis of the fan beam 
produced by the simulator without altering the original 
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f unction of the x-ray source of the simulator, which 
comprises a plurality of individual elements. Each 
individual detector element produces an output signal 
having an amplitude proportional to the energy intensity 
5 of an x-ray incident thereon, the intensity of an 

incident x-ray being proportional to the density of the 
target through which the incident x-ray passes before 
striking the individual detector element. Also provided 
is means for collecting the output signals of each 
10 individual detector element at a plurality of incremental 
angles as said detector array and the source of the fan 
beam are rotated around the target object and translating 
the signals into a back projected computer tomographic 
scan of the target object. 



15 



20 



Translation is accomplished by correcting for the 
spread of the x-rays in the fan beam and for the 
differences in the intensify of the x-rays comprising the 
fan beam depending on the position of the individual 
detector element in the detector array relative to the 
central axis of the fan beam. Translation also involves 
scaling the output signal to account both for the 
relative distance from the source of the fan beam to the 
detector array and for the distance from the fan beam 
25 source to the target object. The attachment also 

includes means for correcting the scaled signal for the 
magnification resulting from the spread of the fan beam 
and for transforming the output signal from each detector 
element into an output signal representing the output 
30 signal that would have been produced by each detector 

element had the incident x-ray originated from a parallel 
beam instead of a fan beam. Finally, translation 
involves converting the transformed output signal at each 
incremental angle of the detector array into a gray scale 
35 value for a picture element having a specific set of 

coordinates relative to the coordinates of the detector 
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array and outputting the gray scale value to an 
appropriate display means for displaying the tomographs 



scan. 



Also provided is a method of producing a computed 
tomographic scan of a target object with an x-ray 
simulator comprising the steps of projecting an x-ray 
beam through a target object, detecting the x-ray beam 
incident on a detector positioned on the other side of 
the target from the source of the x-ray beam, defining 
polar and cartesian coordinate systems to describe the 
geometry of the x-ray beam, target object, and detector 
and locating the position of the source of the x-ray 
beam, the center of rotation axis, and the detector on 
the polar and cartesian coordinate systems for each 
incremental angle as the source of the x-ray beam and tte 
detector rotate relative to the center axis of the target 
object. The output signal from each detector is then 
transformed into an output signal representing the output 
signal which would have been produced by that detector 
had the incident x-ray originated from a parallel beam 
source rather than a fan beam x-ray source and the 
transformed output signal is converted at each 
incremental angle of the detector into a gray scale value 
for a picture element having a specific position on the 
cartesian coordinate system. The gray scale values are 
then output to an appropriate display device. 

in accordance with another aspect of the present 
invention, CT scanning capabilities and panoramic imaging 
capabilities are combined into a single apparatus capable 
of operating in two distinct modes, with only minor 
modifications to the hardware. In a first mode of 
operation (the panoramic mode) , the x-ray source and 
sensor are aligned in a plane, and are rotated together 

~-i^> nA t-f desired, the center of 
in a transverse plane, ir aesureu, 
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rotation can also be translated along a predetermined 
curvilinear path, in order to follow the portion of 
anatomy being imaged, for example, the maxillo facial 
region. In a second mode of operation (the CT scanning 
mode) , the detector and x-ray source are aligned in the 
same plane as the plane of rotation. In accordance with 
the present invention, the eccentric movement of the 
center of rotation can be accommodated so that an 
undistorted CT image is produced. 

Figure 1 is a schematic, sectional view of an x-ray 
simulator having an apparatus constructed in accordance 
with the present invention mounted thereto. 

15 Figure 2 is a schematic of the component parts of 

the apparatus of the present invention and the processing 
of the data by those parts. 

Figure 3 is a flow chart of the processing of the 
20 data in accordance with the method of the present 
invention. 

Figure 4 is a schematic representation of the 
geometry of the coordinate system in which the x-ray 
25 source, target, and detector array of the present 
invention are located. 

Figure 5 is a schematic representation of the 
geometry of the back projected image computed in 
30 accordance with the method of the present invention. 

Figure 6 is a graphical representation of an 
exaggeration of a measured center of rotation shift for 
an x-ray simulator rotating 360°. The off-center shift S, 
35 is plotted as a radius from the center versus the angle 
in degrees. 
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Figure 7A is a graphical representation of a method 
for measuring the center of rotation shift as a function 
of an angle j. 

Figure 7B is a schematic representation of one type 
of error in data acquisition and back projection, which 
results in distortion and blurring of the back projected 
image of a phantom, as a result of the superposition of 
translational and rotational movement of the x-ray source 
of a panoramic x-ray device with respect to the phantom. 
For purposes of clarity, only two rays P l (shifted to 
P 1+ 5) and P 2 (shifted to P 2 +5) are shown having been 
shifted by this eccentric rotation. 

Figure 8 is a schematic representation of the 
divergence of a ray emanating from a point source and the 
different magnification factors that result as the ray 
passes through the object. 

Figures 9A and 9B show alternative methods using the 
spacing element Sr which is utilized during back 
projection as either a constant value (Fig. 9A) for the 
horizontal back projection or a variable value (Fig. 9B) 
for the back projection from the detector line. 

Figure 10 is a representation of the rectangular 
pixel grid containing the circle of image reconstruction 
and showing the sum of the back projection geometrical 
factors at each pixel in the 16 x 16 grid. 

Figures llA and 11B illustrate the filters used for 
processing the output signals of the detector elements 
and the effect of each type of filter on the convolved 
projection data, respectively. 



35 
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Figure 12 is another schematic representation of the 
geometry of the coordinate system in which the x-ray 
source, target, and detector array are located in 
accordance with the present invention. 

5 

Figure 13 is another flow chart of the processing of 
the data in accordance with the method of the present 
invention derived from the geometry of Figure 12. 

10 m accordance with a first embodiment of the present 

invention, and referring to Fig. 1, there is shown a 
schematic representation of an x-ray simulator, indicated 
generally at reference numeral 10. X-ray simulator 10 is 
generally comprised of the x-ray tube 12 used for 

15 radiation therapy simulation which produces a fan beam, 

indicated at reference numeral 14, incident upon a target 
object 16 positioned on table 18. X-ray beam 14 is also 
incident upon the detector array (not shown) of the 
apparatus 20 of the present invention. The apparatus 20 

20 is comprised of the linear detector array 26 (see Fig. 
2), the mounting of which is adapted for fitting into, 
mounting next to, or covering the film holder 22 of x-ray 
simulator 10 normal to the central axis of the fan beam 
14 produced by the source 12 without altering the 

25 original function of simulator 10. The film holder 22, 
or cradle, is mounted on top of the image intensifier 
tube 24. The x-ray beam is collimated by using the 
shutters or collimators (not shown) provided by the x-ray 
therapy simulator 10. The beam 14 is preferably shaped 

30 with a width of approximately one centimeter and a spread 
determined by the length of the detector array 26 of the 
apparatus 20 of the present invention. If the projected 
shadow of the target object 16 to be scanned is less than 
the maximum width of the detector array 26, the spread of 

35 the beam 14 is decreased so as to just encompass the 
diameter of the target object 16. 
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material used is ytterbium-gadolinium oxide, a screen 
material commonly used in diagnostic radiology. Other 
scintillator materials may likewise be used to advantage, 
and the screen may be made of different thickness to 
5 improve photon collection efficiency. 

An external sync trigger circuit 28 is provided for 
interfacing the x-ray simulator 10 with an external sync 
circuit 30 to trigger the timing circuitry of the 

10 apparatus 20 of the present invention when the x-ray tube 
12 starts to emit the fan beam 14. External sync circuit 
30 provides the synchronization signal to the 
preprocessor 32 and determines the integration period of 
the detector array 26. The external sync circuit 30 

15 consists of an input (not shown) to provide software 

control of offset and gain calibration reset signals to 
the preprocessor 32, a continuous clock circuit (also not 
shown) to provide the synchronization and integration 
period for the detector array 26 during gain 

20 calibrations, and a variable clock circuit (also not 
shown) which utilizes an EPROM to provide the 
synchronization and integration period for the detector 
array 26 during data acquisition. The preprocessor 32 
receives the analog signal from the detector array 26 

25 and, under control of the external sync circuitry 30, 
digitizes and outputs the data to the microcomputer 
interface circuit 34. 

Computer interface circuit 34 converts the data 
30 received from the preprocessor 32 to a format compatible 
with the input-output (I/O) board 36 of the particular 
microcomputer utilized for processing the data as 
described below. The microcomputer I/O board 36, after 
receiving the data from microcomputer interface circuit 
35 34, inputs the data to the memory of the microprocessor 
38, which is preferably a Motorola 68030 or Intel 80386- 
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based computer. In a presently preferred embodiment, an 
Apple Mac II Plus (tm) with eight megabyte memory, large 
hard disk, 80-140 megabytes, floppy disk and CRT 
displaying a 480 x 512 pixel image (or, alternatively, a 
5 special CRT to display a 512 x 512 image) having a 25 and 
preferably 30 megahertz internal clock is used to 
advantage. In addition, the output of the microcomputer 
I/O board 36 provides the software control signals to the 
external sync circuit 30 for offset and gain 

10 calibrations, the continuous clock and variable clock 
circuits. The microcomputer 38 performs the data 
manipulation routines described below to produce the CT 
image, outputting the processed data to an appropriate ( 
high resolution monitor 40 which provides a gray scale or 

15 color display of the reconstructed CT image, the film 
recorder (matrix camera) 42 which provides a hard copy 
output (film) of the image on the high resolution monitor 
40 and/or the large capacity storage device 44 which 
stores the image data as well as patient information. 

20 A general schematic outline of the data manipulation 

process is set out in Fig. 3. By the term "data', it will 
be understood that reference is made to the processed 
output signals of each individual detector or photodxode 

25 comprising detector array 26, which is proportional to 
the intensity of the x-ray incident thereon. 
Preprocessor 32 integrates the output signal of each 
photodiode and reads out the resultant signal upon 
external command controlled by microcomputer 38. In so 

30 doing, the data is corrected, smoothed, and filtered as 
shown in the step represented by box 46 in Fig. 3. Tins 
' filtering step is required because the individual 

detector elements, or photodiodes, become defective and 
may malfunction either partially or completely. A 

35 correction algorithm is provided which identifies those 
faulty detector elements through comparison with the 
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neighboring detector elements (the narrow separation 
between detector elements utilized allows the 
generalization that the response for adjacent detector 
elements should be relatively close to each other) , and 
5 replaces the bad data with a corrected, interpolated 
value. A smoothing algorithm is provided for extending 
this method by setting a limit, or preset selection 
criterion, on the variation of adjacent individual 
detector elements and interpolating response values for 

10 those detector elements whose variation from the 

neighboring detector elements exceeds that preselected 
criterion. Recognition of when an individual detector 
element malfunctions is determined experimentally by 
exposure to x-rays both in air and through a phantom 

15 target (not shown) and the history of the detector array. 

Characterizing each step in more detail, the data is 
first smoothed in accordance with the formula 

20 P(2) = (P(D + P(3))/2 [l] 

where P(l), P(2) , and P(3) represent the signal outputs 
from each of three successive detector elements in 
detector array 26. The smoothing is extended by 

25 selecting a maximum variation between adjacent detectors 
based upon the detector response history. For example, a 
cutoff limit of ± C.2 or 0.3 times the adjacenc value may . 
be selected so that if the reading from a particular 
detector element varies from its neighbors by more than 

30 0.2 or 0.3 times the output of the neighboring elements, 
that output signal is replaced by the average of the 
output signal from the neighboring elements on either 
side of the bad element. In the event that several 
detectors in a row have malfunctioned, a scaled averaging 

35 method is used to correct the data. The difference 

between the last good output signal at one end of the 



malfunctioning row of detector elements and the first 
good output signal at the other end is divided by the 
number of detector elements malfunctioning and the 
resultant quotient is used as a constant which is 
sequentially and iteratively added to the last good 
output signal to provide a replacement for the next 
output. The correction is added to this value for the 
next reading and so on, up to the first good output 
signal. For example, given the following output signals, 
110, 125, 255, 255, 255, 255, 140, 141, the correction is 
applied as follows: 

C = (140 - 125) /5 = 3 t 2 ] 

in this example, there are five intervals from the last 
good output signal at one end of the malfunctioning 
detector elements and the first good output signal on the 
other end of the malfunctioning detector elements. The 
corrected output data is 110, 125, 128, 131, 134, 137, 
140, 141. 

Other corrections account for variations in x-ray 
output. These variations can be caused by current and/or 
voltage variations. Metering circuits on the x-ray 
simulator 10 are used to monitor these parameters, and 
the output of those circuits is used as an input to 
processor 32. Current variations are corrected in linear 
fashion: if, for instance, the current drops from 5 mA to 
4 mA, the detector elements are corrected as (5/4) x 
P(i). voltage corrections are non-linear. The 
correction factor is formed as a power function of the 
voltage, i.e., V N , where N is a number in the range 1.5 < 
N < 3. N is determined experimentally and is 
idiosyncratic to the system being controlled. Suppose 
the voltage drops from 120V to 110V and, for that range, 
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it has been determined that N = 2. The correction is 
made as follows: 

P(i)=orr = [120/110] 2 x P(i) [3] 

Additional corrections are performed using 
interpolation procedures to optimize accuracy in the 
calculated intensities used in back projection when 
correlating the detector in the detector array 26 that 
intercepts a back projected, or parallel beam, ray 50 
(Fig. 4) which passes closest to the point in the target 
object 16 being reconstructed. The coordinate system in 
the back projection process is defined from i = 0 at the 
first (-x) detector to N D - 1, where N D is the number of 
detectors in the array 26. For example, using a 2048 
detector array 26 with a detector spacing of x = 0.045 
cm, the total length D L of the line of detectors in 
detector array 26 is 

D L = (2048) X 0.045 cm = 92.16 cm. [4] 

The factor 2048 is used because the detector line is 
based on the spacing of the individual detectors and 
there are 2047 spaces between the 2048 individual 
detectors in a presently preferred embodiment of detector 
array 26. It should be understood that detector arrays 
with other numbers of individual detectors (for example . 
1048) may also be used. The center of the detector line 
is taken as zero. The x coordinate along this line then 
varies from -92.16/2 to 92.16/2, and the x coordinate of 
the i* detector is 

X: = -46.08 cm + (i + 0.5) x 0.045 cm, [5] 



e.g., for the 1135 th detector, 
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= -46.08 + (1135 + 0.5) x 0.045 = 5.0175 cm. 



The coordinate on the detector line where the ray 
passes closest to the reconstruction target point is 
designated L* (back projected) . Once this coordinate is 
determined, it is possible to calculate which individual 
detector will receive the ray passing through I*,. Let 



= intCd^p + 46.08) /Sx) + 0.5] 



[63 



where »int[]" indicates the integer part of the number in 
the brackets in equation [6] and the factor 0.5 is an end 
correction such that in the limit where 1^ is . 
approximately equal to -46.08, the formula will return 
15 the integer value of one (1) . 

Round off error in calculating this integer is 
minimized by linear interpolation. The fractional part 
of Nbp is 



N„p " [Nbp] = fbp- 



[7] 



in the simplest round off, i.e., if f* < 0.5, - [H*] - 
If > 0.5, Nbp = [N^] + 1. For example, if H*« 756.528, 
the detector output taken would be Pl(757) . To be more 
exact, interpolation yields 

pi = (i - f^) PKPV> + f » P1 (£ N * 3 + 1} * 

in the numerical example, Pi - (1 - -528) Pl(756) + .528 
Pl(757) . 

The term "smoothing" (steps 46 and/or 72 of Fig. 3) 
refers to the filtering of the data by both smoothing and 
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ripple filtering. Smoothing is accomplished using a set 
of reconstructed values which may be represented as 

M(1,N) 
M(2,N) 
M(3,N) 



M(N,N) [8] 

Here, N = 2 k where k = an integer, i.e., N - 128, 256, 512, 
etc. For four point smoothing for nearest neighbor (j 
denotes the row and i denotes the column) : 

15 

M (j,i) = [M(j-l,i)+M(j+l,i)+M(j,i-l)+M(j,i+l)]/4 [9] 

The top, bottom, and adjacent left and right neighbors are 
used with equal weight. For five point smoothing for 
20 nearest neighbor, 

M(j,i) = [M(j,i)+M(j-l,i)+M(j+l,i)+M(j,i-l)+M(j,i+l)]/5 [10] 

M(j,i) is added in along with the four nearest neighbors. 
25 For weighted neighbor smoothing with the central value 
given the greatest weight, nine point smoothing is used: 

M'(j/i) = [lxM(j-l,i-l) + 2xM(j-l,i) + lxM(j-l,i+l) 

30 + 2xM(j,i-l) + 4xM(j,i) + 2xM(j,i+l) [11] 

+ lxM(j+l,i-l) + 2xM(j+l,i) + lxM(j+l,i+l) ]/16. 

This weighting system can be varied with different weights 
35 applied to different locations depending upon the 
circumstances as known to those skilled in the art. 



10 
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Ripple filters are used to eliminate the wave 
introduced into the reconstruction by employing Fourier 
methods. The first step in application of such a filter is 
to obtain x-ray transmission values for a water (uniform 
target density) phantom. With known reconstructed CT 
values, these true values are compared to observed CT 
values. The ripple filter factor for each picture element 
is given by the true value/ observed value. To achieve 
correction, the value for each picture element is 
multiplied by the filter factor. There are many other 
different types of filters which may be used to advantage 
depending upon the circumstances as known to those skilled 
in the art who have benefit of this disclosure. Such 
filters might include, for instance, edge enhancement 
filters, high frequency suppression, bone enhancement 
filters, soft tissue enhancement filters, and so on. 

The smoothed, corrected data from step 46 is convolved 
at step 47 and then back projected at step 48. Convolution 
38 and back projection 48 are accomplished in accordance 
with an algorithm which is derived as follows. Referring 
to Fig. 4, there is shown the geometry for the fan beam 14, 
the parallel beam 50 and the coordinate system used to 
describe the algorithm of the present invention. In this 
system, the origin is established at the nominal center of 
rotation 56 with the knowledge that the nominal center may 
vary during the rotation. The x-axis 52 is on a line 54 
through the center of rotation 56 so that the +x axis will 
coincide with the left side of the target object 16 and 
perpendicular to the y-axis 58 which passes through the 
center of rotation 56 on a line from the x-ray source 12 to 
the center of the detector array 26 (which is parallel to 
the x-axis 52) , with the +y axis oriented toward the front 
of target object 16. As shown in stick figure in Fig. 4, 
if target object 16 is a patient, the +z axis is oriented 
toward the patient's head; the coordinate system, when 
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defined in this manner, is a right-handed system. The 
x-ray source 12, is assumed to rotate counterclockwise 
about the z axis with the angle of the x-ray defined from 
the +x axis, and the starting point of rotation is assumed 
5 to be 90°, or the x-ray target coinciding with the +y axis. 
In this fashion, images will be reconstructed that show the 
right side of the subject on the left side of the image and 
reflect the viewpoint of looking at the image from the foot 
of the patient toward the head. 

10 

The distance from the x-ray source 12 to the center of 
rotation 56 is D, reference numeral 60. The separation 
between the center of rotation 56 and the plane of the 
detector array 26 is a distance V, reference numeral 62. 

15 The x-ray source 12 and detector array 20 are incrementally 
rotated about the center of rotation 56 through angles T as 
shown at reference numeral 64 in Fig. 5. The increments 
may be equal in size or unequal. When equal increments are 
used, T - 180°/N or T = 360°/N, where N is the number of 

20 increments, and the system may be rotated either through a 
semi -circle or the entire arc of a circle. Rotation 
through the arc defining a semicircle is commonly used in 
dental panoramic x-rays, and the coordinate system defined 
above is equally applicable in that instance. However, 

25 additional detail is also provided below with reference to 
Figs. 12 and 13. In such apparatus, the center of rotation 
56 is not fixed such that both D and V vary at each angular- 
increment. As set out below, the change in the location of 
the center of rotation 56 is measured and, using the method 

30 of the present invention, the back projection is corrected 
for the resulting known error. 

For equal increments, the index, j, will nominally run 
from 0 to 359 or from 0 to 719 depending upon whether the 
35 projections are taken at one degree or 0.5 degree 
intervals. The number of angles used for acquisition of 
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data is limited by the ability of the data collection 
system and the computer system to acquire and process the 
given volume of input data. The increment between angles 
Tj and T j+1 is denoted by *T and is given by 360/N or 720/N. 

The projection angle, Tj , indicated at reference 
numeral 66, is given, in the equal increment case, by 

•p = go 0 + jST where 0 < j < N - 1. C 12 3 
The x-ray source coordinates for each j are given by 

x.(j) = D.cjjcosTj r i3 3 

y.(j) = D.CjJsinTj. C 14 ^ 

Since the origin of coordinates is at the center 56 of the 
target 16, the nominal distance of the x-ray source to the 
origin is 

D.(j) = [x, 2 (j) + y. 2 (j)]*- 



The coordinates, x.(j) and y,(j) must be known for all angles 
T The x-ray source and detector array are assumed to be 
25 rigidly connected and co-moving in the rotation about the 
target 16. As described above, the angle T 0 = 90° defines 
the initial configuration of the system before rotation, 
e.g., at j = 0. At this angle, T 0 , the coordinates of each 
of the detectors in detector array 26 are: 



x d (0,i) = - D L /2 + (i + 0.5) Si C 15 3 

[16] 

y d (0,i) = - v 

35 where Si = D L /N d , D L is the length of the detector line, and 
N„ is the number of detectors in array 26 where 0 < x < N d . 
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Equations [15] and [16] describe the coordinates for non- 
eccentric motion; as described above, when motion is 
eccentric, both x and y can be measured and input or solved 
for iteratively. The Oth detector is considered the left 
5 most detector on the -Y axis, and the spacing is actually 
the distance between the centers of each detector in 
detector array 26. The factor 0.5 in equation [15] serves 
to define the end point of the ray in the center of the 
detector element. For Tj > 90°, the coordinates of each 
10 detector in detector array 26 are given by: 

x d (j,i) = x d (0,i) cos (j«ST) - y d (0,i) sin (jST) [17] 

y d (j,i) = x d (0,i) sin (j$T) + y d (0,i) cos ( j*T) .' [18] 

In Figure 4, the parallel rays 50 are constructed by back 
projection from the detector array 26. The intersections 
of the parallel 50 and fan beam 14 rays define the circle 
of reconstruction 68. 

20 

In defining the circle of reconstruction 68, the next 
step 70 (see Fig. 3) in manipulation of the data is 
accomplished by microprocessor 38 using the following 
steps. The diameter D re of the reconstruction circle 68 
25 along the line 54 is given by the scaling relationship 

D re /D - D L /W. [19]. 

The radius of reconstruction circle 68 is = D re /2 and W = 
30 D + V defines the distance from the x-ray source 12 to the 
center of the detector array 26. The Jacobian of the 
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scaling transformation from the detector array 26 to the 
line 54 along the diameter of reconstruction circle 68 is 

Mj = D/W [2°3 

The scaling factor M, allows transformation of the data 
taken along the detector line to appropriate values along 
the parallel line 54 through the center of rotation 56 and 
lying along a diameter of the reconstruction circle 68. 

The distance between successive back projected 
parallel beams or detector rays 50 at the level of the 
diameter of the reconstruction circle 68 is the same as the 
individual detector aperture projected back to the center 
of rotation as a length 



5R = Si Mi. 



[21] 



The scaling factor may be writteri explicitly for each row 
100 (see Fig. 8) of pixels in reconstruction circle 68. 
The individual detector aperture projects back to the 
center of rotation as a length 

dV = [D/W]5i C 22 J 

An index c is defined such that c = 0 indicates the top row 
100 of picture elements in the reconstruction circle 68 and 
c = N - 1, e.g., c l5 at reference numeral 104 in Fig. 8, 
indicates the bottom row. Similarly, the index b used in 
equation [38] below indicates the columns 102 of pixels m 
the reconstruction circle 68. Here, b = 0 indicates the 
left most column. The magnification M t defined for the row 
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through the center of the reconstruction circle 68 is 
generalized as 

Magp(c) = [D - y(c) ]/w [23] 
where y(c) = D - R« + (c + 0.5) D w . 



D is the height 106 of the picture element of row c and D,,, 
is the width 108 of the picture element of row c. The row 
10 of pixels lying on the diameter of the reconstruction 
circle 68 therefore have width and height given by 

D pI or D py = D rc /N p [24] 

15 where N P J = number of reconstruction pixels; i.e., 256, 512, 
etc. 

To normalize the projected data from the incident 
beam, fan beam projection data is taken for x-ray 
20 transmission through air and through target objects. Let 

P(j,i) = log (WW [25] 

where is the ith detector reading at angle Tj in air and 
25 r pfi!) is the ith detector reading at angle Tj for x-rays 
transmitted through the patient or other target 16. The 
ranges of the integer indices are as defined above. 

Projection data for each angle Tj has been defined as 
30 p(j,i). This function is re-written in terms of the 
spatial variable s and angle <p in the reconstruction circle 
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68 as p(s,«). The one dimensional Fourier transform of 
p(s,<£) at a given angle <p is written as 

[*p](S,«) = Jp(s,*) exp[-27riSs] ds [26] 



-00 



10 in the discrete form, for each angle T Jf the one dimensional 
Fourier transform of the projection data is 



N ""1 

[5p](S,*) = L°P(k,0) exp[-27TiSs(k)] [27] 
k=0 

where s(k) = -R^k+O.S) *s. Recall that Ss = 2R re /N D and 6s 
= the increment between adjacent rays in the spatial domain 
where s is the spatial domain variable, S is the frequency 
domain variable, i = V=Z in this equation and s = 1/S. In 
this application, s is expressed along the line R (65 in 
Fig. 5) in the reconstruction circle 68, and R is displaced 
on angle *, measured counterclockwise from the +x axis, 
integration is taken over the entire real domain where - 



S < °o- 



The one dimensional Fourier transform of the function 
p (s ,p), written only as a function of s for a given *, is 
30 defined in the continuous notation by 



00 



35 C*P](S) » 



p(s) exp[-2triSs] ds. [28] 



It is necessary to perform a two dimensional inverse 
Fourier transform to obtain the linear absorption coeffi- 
cients, M(x,y), describing target optical densities. 
Coordinates (x,y) are used in the Cartesian spatial domain 
containing the reconstruction circle 68 but defined by the 



68 a« p(s,*). The one dimensional Fourier transform of 
p(s,0) at a given angle <p is written as 

[5p](S,«) = |P(S,0) exp[-27TiSs] ds [26] 
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in the discrete form, for each angle T Jf the one dimensional 
Fourier transform of the projection data is 

N "1 

[*p](S,4) -lP(k,0) exp[-27riss(k)] [27] 
k=0 

where s(k) = -R^Ck+0.5) *s. Recall that Ss = 2^/N, and 6s 
= the increment between adjacent rays in the spatial domain 
where s is the spatial domain variable, S is the frequency 
domain variable, i = V=I in this equation and s = 1/S. In 
this application, s is expressed along the line R (65 in 
Fig. 5) in the reconstruction circle 68, and R is displaced 
on angle 0, measured counterclockwise from the +x axxs. 
integration is taken over the entire real domain where - 
< s < «• 

The one dimensional Fourier transform of the function 
P(s,p), written only as a function of s for a given *, xs 
defined in the continuous notation by 
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[5p](S) = 

-co 



p(s) exp[-27riSs] ds. [28] 



It is necessary to perform a two dimensional inverse 
Fourier transform to obtain the linear absorption coeffi- 
cients, M(x,y), describing target «*W 
Coordinates (x,y) are used in the Cartesian spatxal domaxn 
containing the reconstruction circle 68 but def xned by the 



axis along the center of the projected x-ray beam and the 
orthogonal detector array. The equivalent coordinates in 
the frequency domain are (X,Y) . The projection at angle Tj 
for detector array element i, p(j,i), is expressed as a 
function of the new coordinates. Thus, p(j,i) is written 
as p(x,y) in the spatial domain and P(X,Y) in the frequency 
domain. In practice, the one dimensional transformation 
along the line s is equivalent to the two-dimensional 
transformation for points on the same line at the angle 0. 
Thus, the inverse Fourier transform, T x in this notation, 



is 



( J-'P]( X/ y) = | f | r P(X,Y) exp[27ri(xX + yY) ] dXdY [29] 

J J 

Y. m X. m 

where X ±m , Y ±m are the upper and lower bounds where P(X,Y) 
exists. 

In the polar coordinates (S,*), this equation becomes 

[30] 

[*" l P](s,0) = | |P(S,«) exp [27ri(xS cos* + ys sin*) ] | S | dSd* 

J J 

0 0 

The limits in the integrals are rearranged to 

[30] 

[^»P](s,0)= | |P(S,*)exp[27ri(xS cos*+yS sin*) ] | S|dSd* 

J J 

o -s m 

Note the appearance of |S|, the ramp function, as a result 
of applying the Jacobian to transform the coordinate 
systems. The ramp function acts a filter. Here, the 
frequency is determined by the sampling theorem, which 
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states that where a function h(s) is defined over the range 
2Etc (diameter of the reconstruction circle 68) then the 
Fourier transform of h(s) , H(S), is fully described by 
points SS = 1/2R,, apart (the Nyquist frequency spacing) . 
5 conversely, if the range of interest of H(S) is 2S M , then 
h(s) may be sampled at intervals not greater than s = 1/2S M . 
Thus, in equations 30 and 31, S = 0, S x , S 2 ... S M , and S M 
= 1/2SS, the Nyquist frequency. 

10 The modified projection ray P1(S,*) may be directly 

obtained by applying the inverse Fourier transform to the 
result obtained in Equation [25] or by using the 
convolution theorem. The inverse transform is performed m 
the spatial domain. The inverse transform of P(X,Y) 

15 (written symbolically in the frequency domain as P(J,D) is 
just the projection itself p(j,«). 

The inverse Fourier transform of the ramp function is 
taken by forming an even Fourier series of the ramp 

20 function in the frequency domain and using the coefficients 
of the ascending terms as the inverse Fourier transform in 
the spatial domain. The sine term does not appear because 
the ramp filter is an even function in frequency space and 
the contributions from the negative and positive parts of 

25 the frequency spectrum account for the factor, two, 
multiplying the summation. In the present method, the 
first N D terms of the infinite, even, Fourier series are 
used to approximate the convolution. The modified 
projection ray, P1(S,$) is written in terms of the original 

30 coordinates (J, I) as 

N 

P1(J,D = l P(5,ii) h (|i " iiD C32] 

1,-0 

where h([k|) denotes the kth term in the even Fourier 
cosine series; its use has the effect of yielding a 



circular convolution even though the summation is expressed 
monotonically and linearly. Note that h(k) is independent 
of angle T r It is useful to treat the inverse transform of 
the ramp filter in the light as a continuous function. 
This treatment is accomplished by using the integral rather 
than the discrete transform, or converting the summation in 
eguation [32] to its integral analog. h(k) is expressed as 

h(k) = 25s ||s|cos(s7rk$sS)dS [33] 

where 5s = the increment between adjacent rays in the 
spatial domain. The integral in equation [33] is a 
truncated version of the transform of the absolute filter 
function with cutoffs at S± M . Consequently, no frequencies 
outside the ±S M band will be found in the convolved 
projection data. 

This integral is solved to obtain 

h(0) = l/(45s) for k = 0; [34A] 
h(k) = -l/ff 2 k 2 5S) for k odd. [34B] 

When the projection data is derived from parallel beam 
geometry, 5s is constant. For fan beam geometry, 5s varies 
from pixel to pixel in the reconstruction circle. 

Using the explicit form of h(k) , the modified 
projection data is solved for using the inverse Fourier two 
dimensional transform by writing the spatial domain 
convolution symbolically as 



Pl(k,0) = h(k) * p(k,0) 



[35] 



to solve the inside integral of the double integral formed 

by the inverse, two dimensional transform (the * 

tie oonvolution operation,. The conceptual 

leading to equation [35] is used to derive the hybrid 

rlnsLction algorithm used in the method of the present 

invention. Returning to the discrete representation, the 

Fourier inverse transform of the ramp function, |S|, 



written as: 

N D -1 

II 
a=l 



Pl(k,«, - P(0.«/4«. " (l/^s)|(e,*)/(>c-a) 2 [36] 



where if k = odd, « - 0. 2. 4... and if k = even, a = 1,3, 
S..., so |k - «1 - odd and (k - e) * 0. 

If k = odd, 

P 1 (k,*,.«P(k,*,/(4 5 S,,-tl/(^IP^V (k -° ,2l> - t37A1 
If k = even, 

Pi ( k,«) - <»(*.♦)/«« - u^gh:?' 0 "* 1 C37B1 

These two eolations [37A] and C37B] complete the 
calculation of the inner integral in < 35 ^ 

completing the integration of the outer integral is 
accomplished by calculating the integral over *. This 
process is traditionally called back projection. 

Back projection geometry is shown in Fig. 5. The x.y 
increments are given by: 

«x, Sy = t33] 
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where N p = number of pixels in one dimension, i.e. 32, 64, 
128, etc., and = radius of reconstruction circle 68. In 
similar fashion to equations [15] and [16], 

5 x = -Rje + (0.5 + b)Sx [39] 



10 



y = +R„ - (0.5 + C)6y [40] 
and b, c = 0, 1, 2, — N p -1. 

The linear attenuation coefficient for each x,y or b,c is: 



N-l 

M(b,c) = I Pl(j,i 5 (X,y)) STl [41] 
15 i 5 =0 



where N = number of angles and STl is tt/N for both 180 and 
360 degrees of rotation. Here, i 5 (x,y) is determined at 
20 each angle, and j is the number of the parallel ray 50 at 
the projection angle j which lies closest to the point 
(x,y) . The absorption coefficient is calculated from the 
following relations. 

25 The angle <f> of the ray of interest is defined by 

* - tan' 1 (y/x). [42] 

Because <p is used in conjunction with Tj, <p must be defined 
30 in the same fashion as T 3 , or from 0 to 360 degrees. If x 
< 0 and y >0, then 90° < <t> < 180°. If x < 0 and y < 0, 
then 180° < <t> < 270°. If x > 0 and y 0, then 270° < <p £ 
360°. Let 



35 Lg = W r s sin{[Tj - 0]/[D - r 5 cos (Tj - 0)]}, [43] 

where r 5 = (x 2 + y 2 )* [44] 
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and i,<x,y) - INTCL. + DL/2)/(i + °- 5 >] 

As before, »[] is the -integer part of. The factor 0.5 
appears in equation [45] because when L, ' -°<J 2 ' **l x '*> 

5 must be zero. With L, defined, the number of the detector 
corresponding to L, is determined using equations [6] and 
r 4 5] For example, suppose L, was calculated as 15.463 and 
the length of detector array 26, as it is in a presently 
preferred embodiment, was 69.12 cm with a detector aperture 

.0 of 0.045 cm. The number of the detector lying on that 
point is given by: 

i 5 = HJTKC1S.463 + 69.12)/0.045) + 0.5] = 1112. 

5 This calculation indicates that the ray 50 from detector 
' 1112 passed nearest the point <x,y) that is being 
reconstructed . 

The quantity I* is defined differently for the hybrid 
.0 reconstruction algorithm of the present invention than for 
the parallel beam geometry. At each angle, T„ a line is 
passed from the x-ray target position xT,. *q through the 
point x,y and from there to an intersection with the x-ray 
detector line. Thus, L, is defined as the distance t from 
25 the center of the detector line. For each value of x,y, 
there are N back projections. 

in practice, a rectangular grid with Hp x Mp 
dimensions is reconstructed and all values of absorption 

30 coefficients lying outside the circle of reconstruction 
68 are set equal to zero. Using the rectangular grid for 
reconstruction greatly simplifies any smoothing routines 
when they are utilized. Further, interpolation may be us,* 
when a ray 50 does not pass through a point in the back 

35 projections. 
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Ref erring now to Fig. 6, which is a graphical 
representation of an exaggeration of a measured center of 
rotation shift for an x-ray simulator rotating 360° , the 
off-center shift S i (shown at 76c) is plotted as a radius 
5 from the center of rotation 56 versus the angle in degrees. 
When the x-ray beam which should fall on the central 
detector in the detector array 26 does not pass through the 
center of the circle of reconstruction 68, two errors may 
occur. The first error is a linear displacement of the 

10 detector array 26 parallel to the central axis of the fan 
beam 14. The second error is the magnification or 
demagnification of the distance between detector arrays 
(Si) used in the hybrid reconstruction algorithm described 
above. The latter error occurs when the center of rotation 

15 56 is shifted in a direction perpendicular to the line on 
which the linear detector array 26 is positioned (e.g., the 
detector line) . Further, both errors may occur at the same 
time. 

20 Referring to Fig. 7, a simplified method of measuring 

the shift in the center of rotation is illustrated 
graphically. A round rod (represented schematically at 
reference numeral 74 in Fig. 5) which is relatively opague 
to x-rays is placed at the nominal center of rotation shown 

25 by line 78 in Fig. 7A, and data obtained using the method 
of the present invention is shown as curve 84. An 
eccentricity in the center of rotation causing a shift in 
the direction parallel to the central axis of the projected 
fan beam 14 results in projection data for the rod with 

30 center shown in dotted line 80 with projection data 82. 
The curves 82 and 84 shown in Fig. 7A represent the 
reconstructed projection data taken over each angle T } . 
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From this data, the shift in the number of detector widths 
for each angle is measured. This shift is used to correct 
equation [45] as follows: 



i s (X,y) = INT[ (Lg + D L /2]/i - Bj + 0.5] 



[46] 



Where s- is equal to the shift in number of detector widths. 
However, the maximum shift in either direction must be 
limited so that the shifted detector array 50 at the edge 
10 of detector array 26 still passes through the target object 
16. As shown in Fig. 7B, when the central ray (labeled P 2 ) 
of the detector line does not pass through the center 56 of 
the circle of reconstruction 68, an error is introduced 
which blurs or distorts the image of phantom 78 as set out 
15 above. This error results, for instance, from off-center 
rotation of source 12 and detector array 26 through arc 66 
and 27 (referring to the center point of array 26), 
respectively, or, in the case of a dental panoramic x-ray 
unit, from the superposition of rotational and 
translational rotation, such that the nominal center of 
rotation 56 is shifted to the point 57. The shift xs 
measured and used to correct the baclcpro j ection or 
calculated interatively by subjective evaluation of each 
reconstructed image. For example, if ray P. at angle j is 
25 detected as ray P, + 5, it is designated P, + 5 and used to 
re-set the projection value at P + 5 and so on for P 2 , (P 2 
+ 5), etc. in the application described above, the • 
magnification or demagnif ication of SL is ignored. 

30 In a presently preferred embodiment, however, the 

variable demagnif ication actually present at each row of 
pixels 100 (see Fig. 8) in the reconstruction circle 68 xs 
accounted for in the reconstruction calculation. The 
actual variation of the Jacobian of the transformation xs 

35 used in the method to apply the appropriate value of 
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demagnification, Magp(c), to insert the exact value of the 
back projected data in the calculation. 

X-ray magnification is a harmonic function with the 
mean value M, which can be calculated throughout the 
diameter of reconstruction circle 69 by the following: 



c=N p -l 

m = i/{(i/N.) E[l/Magp(c)]} ; [47] 
c=0 



Fig. 4 shows the calculation of Magp(c) for 16 pixels for 
a circle of reconstruction with U = 100, V - 58, = 21.87 
cm and D„, = 2.73 cm. For example, in pixel (7,9), the 

15 optical density 79.2 (see Fig. 11) is represented as the 
result of the application of projection through thirty 
angles into 64 detectors. Magp(c) is the scaling factor 
indicating demagnif ication caused by beam spread as one 
moves from the detector line towards the source 12. 

20 Magp(c) is measured at the original pixel line, c, with the 
source at 90° or j = 0. 

Magp(c) is the Jacobian of the transformation for a 
coordinate system the origin of which is in the detector 
25 line with points (X',Y') to a coordinate system located at 
the pixel row c with points (X,Y) . For the Jacobian, 



J[(X,Y)/(X',Y')3 



SX/SX' 6X/6Y' 
SY/SX' SY/SY' 



[48] 



J = I Magp(c) 0 I = Magp(c). [49] 
I 0 1 | 

Thus, at the pixel line c, 

X' d(I) = X d(i) Magp(c) C 50 } 



40 
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where i is the index for the detector line end 0 < i . < »o - 
1 starting at the left side of the detector line when j - 
; X correction for the magnification -suiting fro. 
the spread of the x-ray fan beam 14 can he better 

, visualized by referring to Figs. 1, 8, and 9 As best 
shown in Figs. 1 and 8, the spread of x-ray beam 14 is such 
STthe ilge of a target object 16 positioned closer t 
the detector line will be waller than the «age of a 
target positioned farther from the detector line 

0 conversely, in the back projection. 

applied to the data detected at the detector line 

fnseguently, the width of a detector element n. when back 
projected has a dimension Si- = 61 * Magp(=) . 
factor occurs in treating the spacing element ir reference 

. numeral 110 in Fig. 9A, or 112 in Fig. 9B This spacing 
element, as utilized during back projection, is either 
element, latter being 

constant (Fig. 9B) or variable (Fig. 9A) , t* 
preferred because it results in an orthogonal grid of 
pixels as shown in Fig. 8. 

! ° it is Known in the art that convolution filters use a 

ouantity is, which is the spacing between rays, that was 
originally defined for parallel bean geometry. However, in 
^ method of the present invention, the variable spacing 

J5 of the rays resulting from fan beam geometry is accounted 
Without initially using the guantity *r as is used in 
parallel beam geometry. Instead, this quantity is. 
bunted for, in conjunction with HagpCO , in scaling each 
row of pixels as set out above. 

30 The ray spacing appropriate for each pixel is 

accounted for in the convolution using combinations o* four 
filters, ramp, sine. Cosine and Hamming (the functions of 
Lse four filters are shown in Fig. 11A) . The method for 

35 Ting any of these four filters follows that set out above 

^h-^t that is, the filter function is 
for the ramp filter, tna-c is, ux 
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convolved with the function data to produce convolved 
projection profiles at each pixel, and is not repeated 
here. The convolution filters are the Fourier, cosine 
series coefficients and are multiplied by the projection 
data to obtain the convolution profiles (see Fig. 11B) . 
Additional description of the use of such filters for 
similar purposes can be found in such references as S. Webb 
(Ed.), The Physics of Medical Imaging, Philadelphia: Adam 
Hilger imprint by IOP Publishing Ltd. (1988), CHap. 4. 
Each of these filters has certain desired effects, and 
their use is selected to enhance, for instance, bone-soft 
tissue contrast, soft tissue, and so on as described above. 

The geometry of Fig. 7B occurs, in accordance with the 
present invention, in a panoramic x-ray imaging apparatus 
in which the x-ray source and linear sensor have been 
adjusted to be placed within the plane of rotation. 

An exemplary panoramic x-ray image apparatus is 
disclosed, for example, in U.S. Patent No. 5,018,177. That 
apparatus produces digital panoramic x-ray images of a 
target object by impinging an x-ray beam upon an object to 
be imaged and then detecting x-rays that pass through the 
object with a linear detector that produces a pixelized 
electrical signal. The x-ray source and detector are 
aligned with each other, and are rotated in a transverse 
plane about the object being imaged with the center of- 
rotation following a curved path around the object being 
imaged. As the center of rotation traverses the curved 
path, the pixelized electrical signal is integrated during 
predetermined time intervals, and the integrated signal is 
electronically processed to produce a panoramic x-ray image 
of the object. To convert that apparatus from a panoramic 
imaging machine into a CT scanning machine, the linear 
radiation detector, which in this embodiment has 1024 
individual detectors, is aligned horizontally, and the x- 
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ray beam orientation is changed from vertical to 
horizontal. The resultant apparatus, in accordance with 
the present invention, can be easily configured to 
facilitate changing from a vertical orientation for the x- 
ray source and detector for panoramic imaging, or to a 
horizontal orientation for CT scanning. With such an 
apparatus, operation in a first mode produces panoramic 
images, and operation in a second mode produces CT images 
which have been corrected for the eccentrically travelling 
center of rotation of the sensor and detector. 



Fig. 12 presents the eccentric center of rotation 
geometry of Fig. 7B in more detail. Referring to Fig. 13, 
circle of reconstruction 101 has an origin 0 at coordinates 

15 0, 0. X-ray source T is located at coordinates xt, yt and 
x-ray detector 102 begins at point A with coordinates xa, 
ya and extends to point B with coordinates xb, yb. Points 
Kl and K2 having respective coordinates of xl, yl and x2, 
y2, are points lying on circle of reconstruction 101, and 

20 determine a line K that is parallel to detector 102. Point 
D, a point within detector 102, is located at coordinates 
xd, yd. 

Several distances appearing in the geometry of Fig. 12 
25 are also calculated including: dt which is the distance 
between point D on detector 102 and x-ray source T, uc[j] 
which is the distance between the origin O and x-ray source • 
T; ut[j] which is equal to the distance between x-ray 
source T and point U in Fig. 12; ry[k] which is equal to 
30 the distance between the origin 0 and point K0 on the line 
between points Kl and K2; and rx[k] which is equal to the 
distance between point Kl and point K0. It should be 
emphasized that as x-ray source T and x-ray detector 102 
rotate together, with a center of rotation following a 
35 predetermined curvilinear path, the transmission data 
acquired by detector 102 is read out once every 
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approximately 1° of rotation. In this exemplary 
embodiment, an entire rotation spans 231% with each degree 
resulting in a different value for the geometric quantities 
shown schematically in Fig. 12. In addition, for each 
degree of rotation, a number of lines K are calculated, 
each line K being parallel to detector 102. 

Referring now to Fig. 13 as well as Fig. 12, the 
method of the present invention used to correct for gross 
intentional eccentricities in the center of rotation is 
presented in flow chart form. 



Beginning in block 103, for each degree j, the 
detector is read. Then, in block 104, the detector data is 
transformed to produce raw projection data. Control then 
passes to block 106 where the raw projection data is 
convolved with a ramp filter to produce filtered projection 
data. Control then passes to block 107 where coordinates 
of x-ray source T and detector D are calculated using the 
projection angle j and a formula used to correct for the 
shape of the path of the center of rotation. Specifically, 
for each angle j , constant "at" is defined by the equation 
at = (j-26) x 7T/130+7T. Constants cj and sj are then 
calculated as the cosine and sine, respectively, of at. 
25 These quantities are then used as correction factors in the 
remaining calculations. This correction is presented in 
more detail in the source code listing appearing in the 
Appendix. 



Control then passes to block 108 which begins 
calculations involving pixel line Kl to K2. In block 109, 
the magnification along the pixel line Kl to K2 is 
calculated by dividing the distance TK by the distance TD. 
Then, in block 111, the coordinates of the intersection 
points Kl and K2 of line K with the circle of 
reconstruction are calculated. Control then passes to 



block 112 where calculations are performed to determine the 
detector intersection point Dl and D2 from the projection 
rays from x-ray source T through Kl and from source T 
through point K2. Then, in block 113, the detector numbers 
A and B are determined which correspond to positions Dl and 
D2 on detector 102. Then, the functions of nested blocks 
114 116, 117 and 118 are performed whereby the 
contributions of the individual detector elements are 
calculated. Control then passes back to block 108 where 
processing of another pixel line is begun. After all pxxel 
lines within the circle of reconstruction have been 
processed, control passes back to block 103 where the angle 
j is incremented, new detector data is taken, and the 
calculations begins again. After calculations for each 
angle j have been completed, control passes to block 119 
where the CT numbers are calculated. Finally, in block 
121, the results are displayed. 

Fig 13 is a flow chart description of the source code 
of the computer program in accordance with the present 
invention . 

An alternative method used to accommodate eccentric 
centers of rotation will now be described. Initially, 
eccentricities in the rotation of x-ray tube 12 must be 
determined and the individual detectors comprising detector 
array 26 calibrated. A kit (not shown) comprised of 
several rods in a holder, aluminum filter, and a CT phantom 
can be provided for that purpose. Placing the center of 
the phantom for determining motion eccentricities at the 
nominal center of rotation 56, a scan is performed to 
identify any shift in the center of rotation 56. That scan 
also enables verification of the angle encoder output 
whether that angle encoder (not shown) is integrated with 
the apparatus hardware (see Fig. 2) or output from those x- 
ray simulators which include such an encoder. 
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The aluminum filters are then used to determine the 
effective energies and x-ray output for the nominal 
operating voltage of the x-ray simulator. That data is 
used for correcting individual detector output at different 
voltages as at step 46. With a second phantom comprised of 
different materials, scans are performed at different 
voltages. The results are used to set up a Hounsfield 
scale for the different energies. 

Once the initial measurements have been accomplished 
to load the correction algorithms into microcomputer 38 and 
the detectors of detector array 26 calibrated, the 
patient/target object 16 is positioned and the desired 
current and voltage parameters selected. The computer 
program is then initiated to acquire CT data, and on 
computer ready signal, x-ray output is initialized first 
and the x-ray rotation started. For j=0 to N, where j is 
the angle in degrees and N is the number of measurements 
(360, 720, etc.), the computer is signaled to acquire data 
at approximately one degree intervals (or any other odd or 
even increments as selected by the operator) . For i=0 to 
N D -1, when N D is the number of detectors comprising detector 
array 26 and is operator selectable, PP(j,i) is acquired 
where PP(j,i) represents the output of the analog to 
digital converter (ADC) in microcomputer I/O board 36 for 
the detector i at the angle j. Note that the number of 
increments (angles) multiplied by the number of detectors, 
will be the number of individual data points and that, for 
instance, where 0 £ PP(j,i) < 4095, 4095 corresponds to a 
12 bit ADC, 8191 corresponds to a 13 bit ADC, 16,383 
corresponds to a 14 bit ADC, and so on. 

Each of the quantities PP(j,i) is then multiplied by 
a correction number obtained from the calibration procedure 
(step 46) to adjust the data to the values that would have 
been read if the individual detector were operating 
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The aluminum filters are then used to determine the 
effective energies and x-ray output for the nominal 
operating voltage of the x-ray simulator. That data is 
used for correcting individual detector output at different 
voltages as at step 46. With a second phantom comprised of 
different materials, scans are performed at different 
voltages. The results are used to set up a Hounsfield 
scale for the different energies. 

Once the initial measurements have been accomplished 
to load the correction algorithms into microcomputer 38 and 
the detectors of detector array 26 calibrated, the 
patient/ target object 16 is positioned and the desired 
current and voltage parameters selected. The computer 
program is then initiated to acquire CT data, and on 
computer ready signal, x-ray output is initialized first 
and the x-ray rotation started. For j=0 to N, where j is 
the angle in degrees and N is the number of measurements 
(360, 720, etc.), the computer is signaled to acguire data 
at approximately one degree intervals (or any other odd or 
even increments as selected by the operator) . For i=0 to 
N D -1, when N D is the number of detectors comprising detector 
array 26 and is operator selectable, PP(j,i) is acquired 
where PP(j,i) represents the output of the analog to 
digital converter (ADC) in microcomputer I/O board 36 for 
the detector i at the angle j. Note that the number of 
increments (angles) multiplied by the number of detectors, 
will be the number of individual data points and that, for 
instance, where 0 < PP(j,i) < 4095, 4095 corresponds to a 
12 bit ADC, 8191 corresponds to a 13 bit ADC, 16,383 
corresponds to a 14 bit ADC, and so on. 

Each of the quantities PP(j,i) is then multiplied by 
a correction number obtained from the calibration procedure 
(step 46) to adjust the data to the values that would have 
been read if the individual detector were operating 



Thi s correction also remedies the beam 

'"f artifact as discussed above. A reference 

hardening artifact: ^ . . . •_ the 

detector out is also used to correct for ? ™ » «- 
output of the x-ray unit as a function of the angle 3 . 

PCI i) is the input to the hybrid reconstruction 
^/convolution routing 7 *«^^Z 
convolution Fourier art fact filter^ + ^ _ ^ 
described above. For c + o to oy 

tt e bacK projection is done (step 38, and then repeated at 
the next b, next c. 

Reconstructed linear attenuation coefficients are then 
convert to gray scale values using the Bounsfield scale 
1^ fro, scanning the phantoms described above rn 
connection with the setup procedures. A water phantom is 
canned and the reconstructed coefficients are ass^n ed 
„00, with air being 0. Other values, such as for bone 
are then automatically given a value in this scale. In 
Zs fashion, the Hounsfield scale is potentially wider 
rhan that of prior Known C T scales, enabling detect^ 
variations in CT values not currently available By way of 
example, if the scan of a water phantom gives av^ 
0 18 per cm for an average linear attenuation =°« ££lol «* 
f water, 1000 is assigned as the value for water ^d 
1000/ .18 = 555.55 is obtained as the scale factor for the 
unit being used. As is known in the art, 0.36 x 5555.55 
T 2000 is assigned a gray scale of 0 to represent the 
Lsest scale for this system, or twice the density of 
^er and 2 000/ 2 55 - 7.843 is the distribution of the gray 
scale number vs. CT numbers: 

0 CT scale = air = 255 gray scale 
1000 CT scale = water = 127 gray scale 



2000 CT scale - twice density of water = 0 gray scale 

The CT values may be displayed along with gray scale 
numbers is desired. 
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1 in an x-ray apparatus including an x-ray source and a 
substantially linear x-ray detector, a method of operating 
the apparatus in two modes to produce both panorama and 
computed tomographic images of an object, comprising: 

in a first mode, 

positioning said linear x-ray detector in a first 
plane, 

rotating said x-ray detector and said x-ray 
source about an object and about a center of 
rotation in a plane of rotation transverse 
to said first plane while moving said center 
of rotation along a predetermined path, 

electronically processing signals produced by 
said linear detector to produce a panoramic 
x-ray image of said object, 

and in a second mode, 

positioning said linear x-ray detector in said 
plane of rotation, 

rotating said sensor and said x-ray source within 
said plane of rotation about an object and 
about a center of rotation, and 

electronically processing signals produced by 
said linear detector to produce a computed 
tomographic image of said object. 
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2. The method of claim 
second mode, moving said 
predetermined path. 



1, further comprising, in said 
center of rotation along said 



3. The method of claim 1, said detector producing 
pixelized signals proportional to x-rays incident upon said 
detector, said step of electronically processing during 
said first mode comprising: 

periodically integrating said pixelized signals during 
a series of time intervals as said center of 
rotation moves along said predetermined path. 



4. The method of claim 3, wherein said time intervals 
a function of an elapsed time required for said center 
rotation to move along said predetermined path. 



5. The method of claim 1, wherein said step of 
electronically processing during said second mode 
comprises: 

scaling signals produced by said detector in 
accordance with a signal that would have been 
produced by said detector without said object 
being present; 

correcting the scaled signals for magnification 
resulting from spread of a fan beam produced by 
said x-ray source; 

convolving the corrected scaled signals to produce a 
value representing a signal that would have been 
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produced by the detector had the fan beam been a 
collimated beam; and 

repeating said scaling, correcting and convolving 
steps at a plurality of positions as said sensor 
and x-ray source rotate about said object to 
produce a back projected image of said object. 
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6. The method of claim 2, wherein said step of 
electronically processing during said second mode comprises 
correcting the signals produced by the detector for 
distortion from movement of said center of rotation along 
said predetermined path. ■ 

7. A method of operating an x-ray apparatus to computed 
tomographic images of the maxillo facial area of a patient, 
comprising: 

in a first mode, 

generating an x-ray fan beam in a substantially 
vertical plane, 

rotating said x-ray fan beam in a substantially 
horizontal plane about a maxillo facial 
target area while translating said fan beam 
along a curvilinear path, 

detecting an intensity of x-rays in said 
substantially vertical plane after said x- 
rays have passed through said maxillo facial 
target area to produce a detected intensity, 
and 



-57- 



10 



15 



25 



30 



calculating a panoramic image of a maxillo facial 
area of a patient in said maxillo facial 
target area from said detected intensity, 
and 

in a second mode, 

generating an x-ray fan beam in said 
substantially horizontal plane, 

rotating said x-ray fan beam in said 
substantially horizontal plane about said 
maxillo facial target area while translating 
said fan beam along said curvilinear path, 



detecting an intensity of x-rays in said 
substantially horizontal plane after said x- 
rays have passed through said maxillo facial 
target area, to produce a detected 
20 intensity, and 

correcting said detected intensity for said 
curvilinear path to produce a corrected 
intensity, and 



calculating a computed tomographic image of said 
maxillo facial area of said patient in said 
maxillo facial target area from said 
corrected intensity. 



8. An x-ray imaging apparatus for producing both 
panoramic images and computed tomographic images, 
comprising: 
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an x - r ay source positionable to produce an x-ray 
" fan M in first and second orthogonal 



planes; 



linear 



x -ray detector for detecting x-rays 
produced by said x-ray source at a 
substantially fixed distance from said x-ray 
source and for producing signals ind.cat.ve 
of an intensity of x-rays incident on said 
detector, said detector being positionable 
in said first and second orthogonal planes, 

*eans for rotating said source and detector 
of rotation located in said first plane; 

4^ ranter of rotation along a 
means for moving said center or ro 

4-u said source and 

curvilinear path while saiu 

detector rotate; and 

a signal processor for calculating, 

x -ray intensity, a panoramic image of an 
object located near said center of rotation 
when said x-ray fan beam and said detector 
both lie in said second plane, and for 
calculating a computed tomographic .mage of 
an object located near said center of. 
rotation when said x-ray fan beam and said 
detector both lie in said first plane. 



T he apparatus of claim a, said linear x-ray detector 
comprising individual detecting elements. 



The apparatus of claim 9, said linear x-ray detector 
comprising 1024 individual detecting elements. 
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